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Abstract 

It is proved that discrete shock profiles (DSPs) for the Lax-Friedrichs 
scheme for a system of conservation laws do not necessarily depend 
continuously in BV on their speed. We construct examples of 2 x 2- 
systems for which there are sequences of DSPs with speeds converging 
to a rational number. Due to a resonance phenomenon, the difference 
between the limiting DSP and any DSP in the sequence will contain 
an order-one amount of variation. 

1 Introduction 

Consider a strictly hyperbolic n x n system of conservation laws in one 
space dimension: 

u t + f(u) x = 0. (1) 

For initial data with small total variation, the existence of a unique 
entropy weak solution is well known [5], 0. A closely related ques- 
tion is the stability and convergence of various types of approximate 
solutions. For vanishing viscosity approximations 

u t + f(u) x = su xx , (2) 

uniform BV bounds, stability and convergence as e — > were recently 
established in 0. Assuming that all the eigenvalues of the Jacobian 
matrix Df(u) are strictly positive, similar results are also proved in 
PP for solutions constructed by the semidiscrete (upwind) Godunov 
scheme 

Tt U > {t) + ~Kt, [-f ~ ("j-i(*))] =°> Uj(t)=u(t,jAx). 
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In the present paper we study the case of fully discrete schemes, where 
the derivatives w.r.t. both time and space are replaced by finite differ- 
ences. 

We recall that, for the 2x2 system of isentropic gas dynamics, 
the convergence of Lax-Friedrichs and Godunov approximations was 
proved in |7j, within the framework of compensated compactness. Fur- 
ther results have been obtained for straight line systems, where all the 
Rankine-Hugoniot curves are straight lines. For straight line systems 
of two equations LeVeque and Temple JT] utilized the existence of 
Riemann invariants to prove stability and convergence of the Godunov 
scheme. For n x n-systems in the same class, uniform BV bounds, 
stability and convergence of a relaxation scheme, Godunov and Lax- 
Friedrichs approximations were established in (Hj, 0], 53 1 respectively. 
The analysis relies on the fact that, due to the very particular geom- 
etry, the interaction of waves of the same family does not generate 
additional oscillations. 

A key ingredient in the arguments in |5j and pQ is the local decom- 
position of the approximate solutions in terms of travelling waves. To 
achieve a good control on the new waves produced by interactions of 
waves of a same family, it is essential that the center manifold of travel- 
ling profiles has a certain degree of smoothness. We show in this paper 
that this smoothness is lacking in the case of fully discrete schemes. 
As remarked by Serre |13| , for general hyperbolic systems the discrete 
shock profiles cannot depend continuously on the speed cr, in the BV 
norm. In the present paper we construct an explicit example showing 
how this happens. 

Our basic example is provided by a 2 x 2 system in triangular form 

u t + f(u) x = Q, (3) 
vt + g{u) x = 0. (4) 

The characteristic speeds are and f'(u) and the system is strictly 
hyperbolic provided f'(u) > 0. The Lax-Friedrichs scheme for ©-0 
with Ax = At takes the form 

U n+l,j = -(u n<j+1 + U n ,j-l) ~ 7j(f(Un,j+l) ~ f(u n J-l)), (5) 
Vn+l.J = + Vn,j-l) - -(g{u n>j+1 ) - g(u n>j -!)). (6) 

For the rest of the paper we fix a flux function /(it) which satisfies 
f'(u) > 1/4, say, and also the CFL condition |/'(u)| < 1, for all u e R. 
We also take the other flux function g constant outside a bounded 
interval. 

A discrete shock profile (DSP) with speed A for JSJ-J^Jl is a pair of 
functions 

(l7(a:),V(a:)) = (U w (x),V w (x)) 
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satisfying 



TTI u U(x + l) + U(x-l) f(U(x + l))-f(U(x-l)) 
U{x-\) = , (7) 



v{x x) _ V(x + l) + V(x-l) g(U(x+l))-g(U(x-l)) ^ (g) 

We will give a rational speed A = p/q and a sequence of rational 
perturbations e n for which the difference V'^ — y(*+ e ») contains an 
0(1) amount of variation, uniformly with respect to n. The variation 
occurs far downstream in an interval of the form [— Ce~ 2 , —ce~ ]. 



2 Outline of construction 

As a motivation for the later computations we consider the easier case 
of the heat equation with a point-source. If the source acts continuously 
in time and concentrated along the line x = at, then 

V t " V xx = 8 t ,at (9) 

with a > 0. In this case one finds the travelling wave solution 

v(t,x) = (j>{x - at) (10) 

with 

t(y) = J Q Git,y + <rt)dt = { a a _ 1 e * y y <l> (11) 

Here G(t, x) := e~ x2 / 4t /2y / 7rt is the standard heat kernel. Notice that 
the travelling profile can also be obtained as the value at time t = of 
a solution of |j5Jl defined for t s] — oo, 0]. We also have 

<f>'{y) = G x {t,y + °t)dt = { Q e ~ ay ^ y y >° Q (12) 

Next, we consider the case where the sources are located on a discrete 
set of points P n = (n, an), with n integer (the white circles in Figure 1) 

v t - v xx = S, han . (13) 

We again assume that a > and consider a solution of l|13(l defined 
for t e] — oo, 0]. Its value at time t — is now computed as 

v{0,y) = = ^2G{n,y + an). (14) 

n>l 
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For y — ► oo, it is clear that $(?/) tends to zero exponentially fast, to- 
gether with all derivatives. We wish to understand how the oscillations 
decay for y — > — oo, i.e. far downstream from the shock. For y < 0, the 
sum 111 If can be expressed as an integral 

/•OO 

$(y) = V G(n, y + an) = / G(t, y + at) (l + h[ (t)) dt, (15) 
where 

hi(t) = [t]-t+l/2. 

By induction, we can find a sequence of periodic functions h m such 
that (Figure 2) 

= h m (t + 1) , y h m (t)dt = Q, -jjh m (t) = h m -x(t) . 
Integrating by parts and recalling from H15|) we obtain 

S(t/)=/ G(t,y + <7t)(l + — M*)J* (16) 

1 f°° d m 

= - + (-l) m y o —G(t,y + at)h m (t)dt. (17) 

Since 

h m (t) = J hm-i(s)ds, te[0,l], 
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for some point £ m € [0, 1], by induction we find 

\h m (t)\ < 1 Vto>1, Vie 



The identities 

G(t,x) =t- 1/2 G(l,x/Vt), 



Gt — Gj 



imply 



■5—G(t, x) = t"^ 1 )/ 2 • ^—G{\,x/y/i) , (18) 

^ G (M)=t-(-^.^G(l,,M). (19) 

We recall here some basic estimates for the heat kernel and its deriva- 
tives, which we will use throughout the paper. For k = 0, 1, . . . we 
have 



dx k 

Qk 



G(t,-) 
G(t,-) 



= 0(1) 



= 0(1) 



t (fe+l)/2' 
1 

t(2fe+l)/2' 



(20) 
(21) 



In addition we observe that, as y — > —00, the function i 1— > G(t, y + ut) 
becomes exponentially small together with all its derivatives, outside 
the interval centered at \y\/a with width |2/| 5+1 / 2 , for any S > 0. More 
precisely 



sup 

\t+v/<j\<\y\ 5 + 1 / 2 



0(1) 



as y 



-00, (22) 



for some constant c$ > 0. 
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Letting y — > — oo, for every m > 1 the above estimates imply 



m - 




a 





df 



7 G(t,y + at) 



dt = 0(1) ■ y- m/2 . 



(23) 



Similarly, 



d m+l 



dt" 



T G(t,y + at) 



dt = 0{l)-y- (m+1)/2 . (24) 



Since to > 1 is arbitrary, this shows that the function $' is rapidly 
decreasing as y — ► — oo. In particular, taking to = 2 in l|24|) we obtain 
the integrability of <&', hence a bound on the total variation of <£>. 

Finally, assume that the impulses are located not at the points 
Pn = (n, &n) but at the points with integer coordinates Q n = (n, [en]) 
(the black circles in Figure 1) 



v t - v xx 



n, (711 ) 



(25) 



where [a] denotes the integer part of a real number a. 

Again we consider a solution defined for t € ] — oo, 0] and study 
its profile at the terminal time t = 0. Assuming a > 0, a direct 
computation yields 

«(o, y - 1) = := ^2 G ( n > v + M) ■ 



Because of l|23|) . to determine the asymptotic behavior as y — > — oo, it 
suffices to estimate the difference 

K{y) = - $(y) = - ^ [<3(n, y + on) - G(n, y + {an]) 

n>l 

It is here that, if the speed a is close to a rational, a resonance is 
observed. To see a simple case, let a = 1 + e, with e > small. Then 
we can approximate 



-^(y) ~ - ^2 G x (n, y + an) (an - [an]) 

n>l 
poo 

«-/ G s (t,i/ + trt)(et-[et])dfc. 
«/ o 



(26) 



The functions appearing in the above integration are shown in Figure 3. 
We recall that 



G x (t,y + at)dt 



y + at 



exp 



At 



dt = 
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for every y < 0. Set y e = — e 2 . When y ranges within the interval 

I s = [y £ ,y s /2] = [-s- 2 ,-e- 2 /2}, 
the integral in (|26fl can be of the same order of magnitude as 

/ \G x (t,y £ +at)\ dt > c Q y- 1/2 = c e. 
Jo 

Moreover, each time that y increases by an amount Ay — e^ 1 , the 



1 




en - [[en]] 













t 




y l/2 





Figure 3 



phase of the fractional part [[ey]] — ey goes through a full cycle, hence 
the map 

/>OC 

yi-f / G a .(t,y + <rf)(et-[[et]])dt 
Jo 

oscillates by an amount > C\E. In all, we have approximately l/2e 
cycles within the interval I £ . Hence the total variation of the discrete 
profile v]/( 1 + £ ) on I e can be estimated as 

T.V.{¥ 1+£) ; I e ) > c 2 (27) 

for some constant ci > independent of e. We write here \& = *( 1+E ) 
to emphasize that the profile depends on the speed a = 1 + e. By (|27l) 
it is clear that, as e — > 0+, the functions \I>( 1 + e ) do not form a Cauchy 
sequence and cannot converge in the space BV. 

3 Discrete shock profiles for the system 

Since Q is a decoupled scalar equation, the Lax-Friedrichs scheme (JSJ 
for this component admits discrete shock profiles connecting any two 
states, see Jennings ^01- For given left and right states U-,u + wc 
denote by U x (x) the DSP connecting u_ to u + and moving with speed 

, _ [/(«)] 
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We proceed to construct the second component of a DSP with speed A 
by using the discrete Green kernel for © and the Duhamel principle. 
Consider the Lax-Friedrichs scheme for the equation z t — 0, 

— + ! j~ z n,j-l f0Q \ 
z n+l,j - ^ ■ 

We observe that the discrete Green's function K n k for (|28|l is given 

by 

K nk = [ (s)" ((«-fc)/a) forfc = -n,-n+2,...,n-2,n, _ 
\ otherwise. 

Given any DSP for the first equation, a DSP for © is obtained by 
prescribing vanishing w-data at time n = — oo, and then letting the 
it-terms act as a source in © from n = — oo to n = 0. Consider first 
the difference equation 

Vn+lJ = ^( v n,j+l + v n,j-l) + Ipnj, (30) 

where the sources ip n j are assumed given. By Duhamel's principle 
we have that if vanishing data are given at time step — AT, then the 
solution of (|30|) at time step n > —N is 

n-i 

V n,j = /") ^m,j-kK n -l-m,k- (31) 

m=-W-l fcGZ 

To apply this in our situation we introduce the functions tp, H : R — ► K 

by 

and 

:=~ [9(^ + 1)) -5(^-1))] =-\j X ^{s)ds, 

where U = is a scalar DSP for the first equation. In this case © 
may be written in the form l|3Ufl with the sources ip n j given by 

From now on we make the assumption that g(u) is such that ip and 
hence also H have compact support. 

Proposition 3.1. The pair of functions {U^\V^) where is a 
DSP for J3|) and is defined by 

OO 

n=l fcGZ 

is a DSP for the system JSp-©. 
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< Remark 3.1: The proof of this proposition is immediate once it is 
verified that the double sum converges. Since H has compact support 
the convergence follows. > 

We next give a useful integral representation of V^ x \x). For a fixed £ 
we define the function v {x \-\ £) : R — » M by 

oo 

u (A) (x;£) := ^ (JT n _ 1; [ x+A „_ € j + if„_ li j a;+A „_|] +1 ) . 

n=l 

Proposition 3.2. The function V^ x \x) is given by 
V (x) {x) 

Proof. 

oo 
n=l ' 

oo 

= -^EE 



V(0» (A) (*;£K- 



(32) 



/•a; — fc+An+1 



/•re — fc+An 



n=l k£Z _ 

oo 



71=1 keZ 



' x — k-\-Xn 
px— /c+An+1 

'x — fc+An 



-fc+Ara-1 



/■x— fc+An 



-fe+An-1 



^ 00 " f>00 poo 

-~E / ^(0^ n 

„=1 L^-oo J-oo 
^ />oo 00 

/ <K0 E (-Kn-l.ts+An-a+l + ^n-l,[a;+An-£]) 
Z J -°° n=l 

-t /inn 



di 



2.-r 



□ 



4 Approximation of the DSP in terms of 
the heat kernel 

We will compare solutions of the Lax-Friedrichs scheme with certain 
solutions of the heat equation. As a first step we approximate the 
discrete Green's function K n ^ using the heat kernel 

1 _x£ 
G(t,x) = — =e 4t . 
2v^rt 
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We use the following notation (see [S]): 
a k (v) :■■ 



1 \ 2v 



2) \v + k 
By Stirling's formula we have 

a k (v) = hj\f(hk) ■ exp(ei - e 2 ), 
where _ 

/2 1 
-, Af(x) = —j=e 2 

and the errors £i, £2 satisfy 
3fc 2 A; 4 

— — t < £1 < —5-, provided Ifcl < v/3, and e 2 = Oil/v). 

There are two cases to consider depending on whether both n and k 
are even or both are odd. 

• Case 1: n = 2m, fc = 21. In this case we have, 

K n ,k =cn(m) = 2GfQ,fc) -e £1 - £2 . 

• Case 2: n = 2to +1, fc = 2/ + 1. In this case we have, 

K, h k = ——rarfm) 
n + k 

n + fc 



= 2G(-,fc -e £l - £2 -J -e 2 ™ 

V 2 ' / n + fc V n - 1 

We will use this approximation only in the case when |fc| < 0{l)n 1 ^ 2+s 
where < S -C 1. It follows that in either case we have 

e £l - £2 = 1 + 0{l)n- 1+u . 

An easy calculation shows that under the same condition on fc, 



e m-nt^=Tj = 1 + 0(l)n- 1+5 . 



n + fc y n — 1 
Recalling H2()|) . summing up we have the following. 

Proposition 4.1. For n > 1 and /or |fc| < O^n 1 /^ 5 , with < 5 <C 
1, we ftai/e 

n ■ N -3/2+5 



K n>k =2G{-,k) +0(l)n-*' z+d . (33) 
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For a given speed A > and for 6 G]0, l/2[ we define the time 
interval 

I(y;X 7 S)=[l- y ^,l +l //^]. (34) 

We will make repeated use of the fact that, for y < 0, the indices 
outside 7(|y|;A, 5) contribute exponentially little to the sum 

oo 

S (V) '■= E K n,lv+XnJ- 
n=l 

Proposition 4.2. For y < 0, A > and S G]0, l/2[ we have 

E K n<ly+Xnl <0(l)e- c ^ <X ' !) , (35) 

ngJ(|y|;A,«) 

for some positive constants C(8, A) and c(A, 5). 

Proof. For notational convenience assume that A = 1, the case A < 1 
being similar. First of all we have K n ,[\y+\n]] 7^ iff n > We divide 
the above sum in four parts, where n ranges over I\ = 1 2 ' 2^ ] > ^ 2 ~ 

] ^ \y\ - \y\ 1/2+s [,h = ] \y\ + \y\ 1/2+s , 2\y\ [ and h = [%|, +oo [, re- 
spectively. If n G Ii (and there is at most one such n) then [y+ra] = — n 
and n = 0(|y|), so that K n ^ y+n ^ = = 2~ c 'l a l is transcendentally 
small. For the remaining indexes we can use the following estimate 
obtained by Stirling's formula 

n\ Cn n+l / 2 

, < 77X7777 77 tXTTTT- for < fc < n. (36) 

k) ~ k k+1 / 2 (n- k) n - k+1 / 2 



If n G 7 2 , from J2U) and |J3BJ| it follows 



#„,[b+n]] < y „ (2ra _ y)2 n-» = : ™- 

A calculation shows that n i— ► y) is increasing on 7 2 and that 

FxO/) := In fmaxF^j/)) = lnF(|y|-|y|* + «) = _|y|M +0( |j,|3 5 - 1/2)< 
\nei2 J 

The case n £ ^3 is treated in a similar way. It follows that 

E K nAy+nl <0{l)\y\e- c ^\y\ 2 \ 

n£l 2 Ul 3 

Finally let n G -Z4. Since the map y 1— > F(n,y) is increasing when 
|y| < n, then F(n,y) < F(n,n/2) < (2/3) n / 2 for n G Z4. Thus 

2\ ™/ 4 l y l/ 2 

3 



E = E(i) = ^ 



This completes the proof. □ 
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In the following we will also need an analogous result for the heat 
kernel G. 

Proposition 4.3. Let y < 0, A > and S €]0, l/2[. Then, outside the 
interval I(\y\', A, 5) the integral of G as well as any of its derivatives is 
transcendentally small, i.e. for any k > 0, we have 



j 



Qk 

Ox 



k G(t,y + Xt) 



dt < 0(l)e 



-C(6,X)\y\<^ 



/M+\/(|y|;A,<5) 

for some positive constants C(S, A) andc{\,8). Moreover 



E 



Qk 

dx k 



G (n,y + An) 



< 0(l)e 



E 

n<tI(\y\;\,S) 



dxk G(n,ly + Xnj) 



<o(i)e- c ^y^ x ' s \ 



(37) 



(38) 



(39) 



Proof. For simplicity, assume A — 1. Concerning the integral, we have 



v\-\v 



1/2 + 5 



dt = O(l) 



\y\-\y\ 1/2+s 



I 1 2-5 /c 



0(l)-|y|e 



-|y| 2 V5 



/ 



|»|+|i/|V>+* 



eft = 0(1) 



ir + 


rl 




•%l J 



"4(t+|v|) rfr 



< C\y\e-^ 2S / 8 + / e" T / 8 dr = 0(l) e ~ c W. 

Hence (|37|) follows. The other two estimates can be obtained from 
(HSJl and 111. □ 

As a consequence of the previous propositions, we are authorized 
to add or subtract the tails of the integrals/sums, introducing an error 
which is exponentially decreasing with y. This will be frequently and 
tacitly used in the following. 



4.1 Estimates on the approximation 

In order to estimate the variation of the second component of the 
DSPs we will need that v^ x '(x; £) is close, within acceptable errors, to 
the corresponding function defined in terms of the heat kernel. This 
function is given as 

oc 
n=l 
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where 

z n = z n (x,X) :=x + \(n+l)-£. (40) 

We will need to carefully keep track of the dependence of on A. 
We proceed estimate w (A) . Using JTSJ, l|23)l. PD)l and (J2U, for all 

m > 1, a; < and £ in the support of ■0, writing z = x — £ + A, we 
have 



n=l 

= 2 E G 2 £{ G - G (?' w)} 



n=l 

2 



A . o(i)M 



-2 J] •((z„)) + 0(l)n- 3 / 2 ((z„)) 2 } 

ne/(|*|;A,i) 

| -2 ^ ■((*,)) + 0(l)\x\- 1+s , 



(41) 



where ((a)) := a — [a] is the fractional part for any real number a. In 
this calculation we have used that \I(\z\; A, S)\ — 0(l)|x| 1 / 2+<5 and that 
n = 0(l)|i| when n € I(\z\; A, (5). 

The case of rational speed The estimate (|4*T1) is valid for any 
speed A. Now suppose that A is a rational speed, 

A = -, P,geN, 
q 

with p and g relatively prime, and consider the sum on the right-hand 
side of (|41() . By writing n = mq + j with m > 0, j 6 {0, . . . , q — 1}, 
Taylor expanding about the points (t m ,x m ) where 

t m = mq/2, x m = x m (x, A) := x + Xmq - £, 

and using the formula 

£(( 2 + f))=fe)) + ^, (42) 
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from l|2rj)) and we obtain, for z n = z n (x, A) and z = x — £ + A, 



XI G * (^n) (00) 
ne/(|s|;A,4) 



m>0 3=0 



™9+3'e-f(kl;A,<5) 

9-1 



E E { G * + su p(I G **I + \G xx \)0{q)\ ({z mq+J )) 

m>0 3=0 ' ** 

mq+jeI(\z\:\,S) 

y ] Gx ( _ 2 _ '' Z '" 1 ) ^ ^ (( z m<?+j)) 
m>l " 3=0 

+ 0( g )|x| 1 / 2 + 5 (|xr 2 + |xr 3 / 2 )+e- c l a: l (43) 
9-1' 



((g(* -0)) + V E G * (ir> + °(«)l= 



-1+5 



In this calculation we have used that |/(|z|; A,J)| = O(l) 1 / 2+<s and 
that n = 0(l)\x\ when n € /(|z|;A,(5). Analogously to (|24|l . one can 
prove that 

m>l 

for every integer M. We thus have that 

X G x (CO) 

nGJ(|z|;A,(5) 

= {««(* - od + ^y 1 } E G * (i 2 .^) + G (9)i-r 1+5 

^ J m>l 

=0(g M )|a;|-^ 1+5 + 0(<?)|z|- 1+5 . (44) 

From l|41|) and l|44|l we conclude that when the speed A = p/q is a 
rational we have 

u;W(a;0 = \ + 0(q M )\x\- M ^ 1+5 + 0(q)\x\~ 1+5 . (45) 
A 

Due to the dependence on the denominator q, equation (|45|l is not 
useful for computing the variation of differences w^ ' — w 1 -^ as A — > A. 
Instead, (|45|l will be used for a fixed reference speed A, while we need 
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an alternative analysis to estimate w^ x \ where A = A + e is a small 
perturbation of A. 

For this we return to the right hand side of (|43[) . with x rn = x m (x, A) 
and z n — z n (x,X). We establish the following technical result. 

Proposition 4.4. Let A = p/q and A = A + e, with \e\ <C 1. Then for 
z — x — £, where £ lies in the support of tp, there holds 



E G * 



(~Y' Z + ~ Xm i) ■ E(( z + m i £ + ^)) 

G x (— ' i z + Asg) • ((g(z + s<?e))) ds+ 

+ 0(g)|x|- 1+ * + Odelg 2 )^!" 1 /^ (46) 



Proof. We compare both the sum and the integral to the same integral 
and show that in each case the error is 0(q , )|a;| _1+ ' 5 + 0(|e|g 2 )|x| _1 / 2+ ' 5 . 
Let dj = z + Xj and consider the following difference 



E Gx {~Y lZ + ^ mq ) ' E(( z + mq£ + ~ 



m>0 



J = l 



< 



< 



poo Q 

- G x ,z + Xsqj -^((flj -,+sqe))ds 

e/; +1 h(^.-h- 

m>CT m 

9 

-Ga, (y, z + Asg) | 53(( a J + mge)) ds 

/•m+l ,„„ N f 9 

5^ / G x [-f,z + Xsq) I ^fe + mge)) - (( 

5] 53 (sup (|G at | + \G xx \)0(q)\ ((oj + mgej + e" ^!- 

m>0^ — 1 



Oj + sqe)) > ds 



mg+jeI(|z|;A,<5) 

, 0(1) ' 



■r 



E E 

i =1 mge/(|z|;A,<5) 



m+l 



| ((% + mge)) - ((%• + sqe)) | ds. (47) 



The first sum on the right hand side of l|47|) is bounded by 0(<?)|ir| 1+s . 
To estimate the second sum, we divide the interval 7(|z|;A, 5) into 
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r = 0(\s\q-\I(\z\;\,6)\) = 0(\s\q) \x\V 2+s intervals Ji, J r of equal 
length l/(|e|g), and re- write the sum over m as 



Y / | ((Oj + l s il £ )) - (( a j + st l £ )) | ds - 
k=l Jjk 

The integrand is bounded by \e\q except on a sub- interval of length 1 
where it is O(l). Hence, for each k, the integral over is bounded 
by order one. It follows that the second sum on the right hand side of 
(|4*7|) is bounded by 0(|e|g 2 )|a 



-1/2+5 

On the other hand, using (|4^|) . (jT2|) and JSUJ), we have 



G x (y, 2 + ^ s <?) ' {(q(z + s ?£))) ds - 

p oo Q 
poo ( 1 

^ G x (^,z + \ S q)-U q {z + sqe)))-Y,{ 



cij + sqe)) > ds 



i=i 



ds 



J Gx (y' z + As<7 ) ' ^ + Aj + s?e ^ ~ + sqe ^) ds 



O(l) 



?W j=1 J/(| 2 |;2A,6) 



((z + Aj + 2et)) - ((cij + 2et)) dt + 0(l)\x\ 



-M 



0(1) y 



((v))-((v + sj))\d v + 0(l)\x 



-M 



(48) 



where J(\z\;j,e,X,S) = z + Xj + 2el(\z\; 2A, S). Each of the integrals 
in the sum is therefore of order 0(e 2 q)\z\ 1 ^ 2+s . Finally, using that 
z — x — £ and that £ varies in a compact interval, from l|47|) and (|48|) 
we thus obtain □ 
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Using H43|) and l|41|) for speed A together with Proposition ^. 4l vields 
0=* 2£G.(^,* + Amq) • + mqe + 

m > j — 1 

+ 0(g)|a;r 1 + 5 

2 f 00 / s \ 
= ~ - 2 / G x (-|, 2 + AsgJ • ((g(z + sge))) 

+ 0( <Z )|a;r 1 + 5 +0(|e|g 2 )|xr 1 / 2 + 5 

2 2 f 00 /t - \ 

G x [ -,z + Xt) ■ ((qz + qet))dt + E(\x\;e,q), 



A g7 \2 
where z = x — £ and 

We proceed by simplifying this expression. Put y = — x 3> 0, such that 
z + Xt = — y — £ + Ai, and introduce the coordinate r by 

Thus T-values outside an interval [— Cy s , Cy s ] corresponds to i-values 
for which the contribution to the integral is exponentially small. A 
simple calculation now shows that 

G x (* z + Xt) = -^-re^ 2 / 2 + 0(l)y" 3 / 2 +*. 

Substitution into the expression above yields 

w Cx \-y;0 = l+E(y-e,q)- 
A 



1 / 2 f _„-t 2 /2 ft 1 £ (,. _ IV 



\r\<Cy s 



(50) 



Recalling (|4"5|) we get that for a rational speed A = p/q and a pertur- 
bation A = A + e, there holds 

(as;0 - (z; = 1 - £ + e, «,)- 
A A 

^1 /" _,-r 2 /2 /T^ , ? e ^. _ fy 



1/2 / » is ' ~ i y • \ ~ 

qy ' J \ A V V a 

\T\<Cy~ 



(51) 
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where C\ is a positive constant. We will use this estimate for rational 
speeds A converging to a fixed rational speed A. The relevance of the 
estimate is that it does not depend explicitly on the denominator of 
the approximating speeds. 

4.2 Estimating v (A) — 

In this subsection we show that is a sufficiently good approxima- 
tion of v( x \ We have the following estimate which actually holds for 
any real speed A. 

Proposition 4.5. Provided £ is in the compact support oftp, the func- 
tions (x; £) and w^ x \x; £) satisfy 

v^(x-0-wW(x-t)= X) 2G.g,z„)+0(l)M- 1+25 , (52) 

ne/(| 2 |;A,<5) 

when x < 0, where z = x — £ + A and z„ = z.„ (a;, A) is given by 
Proof. Observe that 

G [*J + l) = G (|, [z„]) + G* Q, z„) + 0{l)n-V\ 

By Proposition ^. 2l onlv indices n € /(|^|; A, <5) in the sums defining i/*) 
and vj( x ) are significant. If n S /(|^|; A, 5) then n = 0(l)|x| while [z„] 
is 0(l)|a;| 1 / 2+ ' 5 , so that we are within the validity of the approximation 
(|33|l . Since exactly one of n — n — ([z n ] + 1) is even, we have 

v W (x; = X (*M*J + *».I*J+i) + 0(l)e- c|a;| 
nel(|z|;A,i5) 

= X 2G(|,[z„])+ X 2G(|,[z„] + l 
»el(|z|;A,<5) ~ neJ"(|z|;A,5) 

ti-I» n ]evwi »-IJ„lodd 

+ X 0(l)n _3/2+tf + 0(l)e- c M 

nei"(|z|;A,i5) 

= E 2G(|,[z„])+ X 2G *(I" 
nei"(|z|;A,<5) ~ neJ"(|z|;A,5) 

n-[*„]odd 

+ X 0(l)n- 3 / 2+5 + 0(l)e- c l a; l 
neJ"(|z|;A,5) 



= u;< A > (*;£) + E 2G X (-, z„) + 0(l)|x|- 1+25 , 

ne/(|z|;A,<5) 

7l-[z„]odd 

since J(]^[; A, <5) has length 0{l)\x\ l / 2+s . □ 
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To estimate the sum on the right-hand side of (|52|l we will make 
particular choices for the speeds A and A. In what follows we will fix 

A 4 ~ x = 2fcTT' ^ 

such that e = — 4fc _i 2 ■ In order to satisfy the non-resonance condition 
of Majda and Ralston ^2] we will let fc be an even integer. We have 
the following key estimates that are independent of k (and e). 

Proposition 4.6. With A and A given by l,5,'j|) we have 

G * (f >*«0> A )) = iS. ( 54 ) 



n€/(|z|;A,*) 

n-[z„(x,A)]odd 



and 



E ^g,^,A))=^. (55) 



nG-f(|z|;A,<5) 

n-[z„(x,A)]odd 

Proof. We give the proof for l|55[) . the proof of (|54|l being similar. Given 
£ and x, we set z = X — £ + A such that z„(x, A) = z + An. We denote by 
Af the set of the integers n contributing to the sum H55J1 ■ Let a fc-block 
denote a half-open interval of length 2k + 1 consisting of k consecutive 
subintervals of equal length 1/A = 2 + 1/fc, each on which the function 
B(s) := \z + As] is constant. Without loss of generality we can assume 
that I(\z\; A, S) is exactly partitioned into finitely many /c-blocks. A k- 
block contains 2k + 1 integers and each subinterval has length 2 + 1/fc. 
It follows that in each fc-block there are fc — 1 subintervals containing 
exactly two integers, and one subinterval containing three integers. 
Since B(s) is constant on each of the subintervals, the function n i— > n— 
[z + An] takes on one even and one odd value on each subinterval that 
contains exactly two integers. As fc is even two consecutive fc-blocks 
will contain exactly 2fc + 1 integers n in J\f. Thus, of two consecutive 
fc-blocks, one contains fc, and the other fc + 1, of integers n in Af. We 
observe that in the sequence of these indices n, the elements are at most 
a distance 3 apart from each other. Finally, since fc is an even integer, 
in all subsequent unions of two consecutive fc-blocks, the distribution 
of the indices n € TV is the same. We can therefore define a map n that 
maps Af bijectively onto the regular grid of even integers in I(\z\; A, 5), 
in such a way that \n — fi(n)\ < 3. The sum in (|55|l can therefore be 



19 



estimated as follows, 



J2 G x (J, z n (x,\)) = J2G X (^,z« n) {xA)) + 
= ^ G x Q,2„0,A)J + 

ne2Zn/(|^|;A,i5) 

+ 0(l)sup(|G srt | + |G? xx |).|I(|z|;A,5)| 

_ 0(1) 



where we have used JUJl, and JUJ. □ 

<i Remark 4.1: We note that the complexity of the preceding ar- 
guments is essentially due to the fact that we are working with the 
Lax-Friedrichs scheme. The same computations would be significantly 
simpler for e.g. the upwind scheme, in which case the complication of 
even and odd terms do not occur. t> 

Combining Proposition 14.51 and Proposition 14.61 we conclude that, if 
the velocities A and A are given by J53J), then 

u« (x; £) - v& (x; = W ™ (x; £ - w<*> (x; + . (56) 



5 Variation of - 

In the remaining part of the paper the velocities A and A = A + e are 
given by (|53|l . Recalling that ip has compact support we conclude from 
(E3), (|2)J) and JSU) that 

V™ (x) - V™ (x)=--~ ^(0 dti + E(y, e)+ 
AA J — oo 

oo 

-oo \r\<Cy s 

= A{e)+E{y,e) + -%( Te~ T '' </2 h(r;y,s) dr, (57) 

V 1 J\T\<Cy* 
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where y — —x and where we have introduced 

poo 

E(y,e):=E(y,s,2) 



A{e) x=-X /"%(£) d£, 
AA J — oo 



h(r;y,e) 



£>(* + *-tH/D)* 



Making the change of variables 77 = 2£ and denoting 0(77) = ip(rj/2), 
13 = 2A/A, and 7 = 2/A 3 / 2 , we have that 

1 f°° 

h{r; y, s) = - y 0(77) ((77 + /9?7 + ■ye^/yr)) dn. (58) 

We proceed to use Ij57(l to show that the function — contains 
an O(l) amount of variation on an interval of the form 

J(e) := [-Ce- 2 ^ 1+2S \-ce- 2 ^ 1+2 % 

Recalling g§) we see that E(y, e) is of order 0( e 2 ( 1 -f)/( 1 + 2 ' 5 )) on J(e). 

As a first step we consider the limiting case where <j>(j])/2 is the 
Dirac delta function centered at 77 — 0. In this case 

h{T;y,e) = /i Q (r;j/,e) := ((/% + 7 E \/y r )) > 

such that 

vW(a:) - ^(z) - A(e) + E(y, e) + ^H {y;e), (59) 

where 

H {y;s):= [ re^ 2 {{0y + -ye^yr)) dr. (60) 

J\r\<Cy s 

Lemma 5.1. There exist an 0(e~ 1 ^ 1+2S ^) number of points y\, ...,yL 
in —J{e) with 

f -0(1) /or ?i odd, , 

H (y n ; e) = i 2 s/(i+26)\ r 61 
L O[e M,( ^ >) for n even. 

Proof. Set 

Vi := [ £ - 2/(1+25) ]//5, 



and let 

Vn ■= VI + \ { f»[e- 1/(1+M) l " " ' 
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for n = 2, . . . , [e 1 /( 1 + 25 )]. These choices imply that 



fyn + l£y/lM T )) | T=0 = ((/%")) 



for n odd, 
1/2 for n even. 



We observe that an O(l) change in the constant C in H6Q(I induces an 
exponentially small error (with respect to y). In order to simplify the 
computations we make the following choices for the constant C in l|(jU|) , 



s+Tr> = 0(1) for n odd, 

- — jqrrra = 0(1) for n even. 



Thus, up to an exponentially small error, we have for n odd 
H (y n ; e) = / re _r2/2 ((/3y n + je^/y^r)) dr 

J\r\<C n y & n 
rO 

re- T ^ 2 = -0(1), (62) 

— OO 

while for n even a similar argument yields 

H Q (y n ;e) = 0(l)e 2S /^K 

□ 

For a fixed, small e it now follows from Lemma 15.11 and H59J1 that 
the function (x) — (x) alternates between values 

A(e) + (l) e 2(l-5)/(l+25) _ (i) e i/(i+a*) ) 
for x — —y n , n odd, and 

A(e) + 0(l)e 2 ^^ 1+2 ^ +0(l)e, 

for x — —y n , n even. Provided e is small enough, this implies that in 
every interval [— y n +i , — Un] , the function V^ x \x) — V^ (x) contains an 
0(l)e 1 ^ 1+2<5 ' ) amount of variation. Since there are an 0(l)e _1 ^ 1+2(5 - ) 
number of such intervals, this argument shows that, in the limiting case 
where 4>{rf)/2 in (|58|) is a Dirac delta, the function V^(x) — V^'(x) 
contains at least an O(l) amount of variation on J(e). 

It remains to argue that the same result holds whenever <p(rj)/2 
is close to a Dirac delta function. For this it is sufficient to show 
that the difference \H(y;e) — Ho(y;e)\ can be made arbitrarily small 
independently of the 0(l)-estimate in Ij62[) . This will be accomplished 
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by imposing that the support of tfi (hence ip) be sufficiently small. We 
have 



\H(y;e)-H (y;e)\ 



< 



r\<Cy 6 



K 



(T){((F(v;e,T)-r,))-((Ftoe,T)))} 



dr dr], 



where F(y; e, r) = Py + je^yyr. Since e^fy is O(l) it follows that there 
is a constant C such that 



(T){({F(y;e,T)- V ))-({F(y;e,r)))} 



K 

\r\<Cy s ' 

We thus have that 

"+QO 



dr < Cr). 



<t>(v)v 

-oo 



dr), 



which can be made arbitrarily small, compared to the 0(l)-estimate 
in (fr!2*|) . by choosing <f) sufficiently close to a Dirac delta. 
We have thus proved the following theorem. 

Theorem 5.2. Discrete shock profiles for the Lax-Friedrichs scheme 
for strictly hyperbolic systems of two conservation laws of the form 

do not depend continuously in BV on their speeds. More precisely, 
one can find a sequence of rational speeds X n converging to X 6 Q, for 
which there are discrete shock profiles ^ n and ^ of speeds X n and X, 
respectively, and such that 

T. V.(tf n - *) = 0(1), 

independently of n. 

< Remark 5.1: From (O, lf5I|) and we observe that with 
our particular choices of speeds we have 

^)(, + Ax)-yW(,) = ^, 

for a; < and Ax = 0(1). It follows that an 0(1) translation of V w 
relative to changes the total variation of their difference only by 
0(e) in the region J(e). > 
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